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\ Abstract. We prove a large deviation principle of Freidlin-Wcntzell's type for the multi- 
valued stochastic differential equations with monotone drifts, which in particular contains 

' a class of SDEs with reflection in a convex domain. 

CIh' 

1. Introduction 

Consider the following multivalued stochastic differential equation (MSDE in short): 
dX(t) G b{X{t))At + a{X{t))<lW{t) - A{X{t))<lt 



(N 
(N 



X(0) = xe D{A), 
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where A is a multivalued maximal monotone operator, which will be described below, 
W{t) = {W''(t),t ^ 0, G N} is a sequence of independent standard Brownian motions 
on a fihered probability space {Q, P; (J^Ot^o), 6 : M"^ ^ M"" and a : M"' ^ M"" x are 
^ ■ two continuous functions, stands for the Hilbert space of square summable sequences 
of real numbers. 

This type of MSDE was first studied by Cepa in [SI E]. He proved that if b and a are 
^ ' Lipschitz continuous, then there exists a unique pair of processes {X{t),K{t)) such that 
b '. 

CZ ft ft 

X{t) = x+ / b{X{s))ds+ / a{X{s))dW{s) - K{t), 



where K{t) is a finite variation process (see Definition 12.31 below for more details). Re- 
cently, Zhang [14J extended Cepa's result to the infinite dimensional case, and relaxed the 
Lipschitz assumption on b to the monotone case. It should be noted that when A is the 
subdifferential of the indicator function of a convex subset of R'", the above MSDE is the 
same as the usual SDE with reflecting boundary in a convex domain (cf. [HE]). Moreover, 
since the subdifferential of any lower semicontinuous convex function is a maximal mono- 
tone operator, Cepa's result can also be used to deal with the SDE with discontinuous 
coefficients. It is well known that there are many literatures to investigate the SDEs with 
reflecting boundary since the solutions of a large class of PDEs with Neumann boundary 
and mixed boundary conditions can be represented by the solution of such SDEs (cf. p^). 



Key words and phrases. Multivalued stochastic differential equation, Maximal monotone operator, 
Large deviation principle. 
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We now consider the following small perturbation of Eq.([T]): 



r dX'{t) e b{ X%t))d t + ^a{X'it))dWit) - A{X%t))dt, 

{ X'{0)=xe D{A), e G (0,1]. ^ ^ 

The solution of this equation is denoted by {X''{t,x), K''(t,x)). We want to establish the 
large deviation principle of the law of X^(t,a;) in the space S := C([0,T] x D(A); D{A)), 
namely, the asymptotic estimates of probabilities P{X'^ e F), where F G i3(S). 

In [I], Anderson and Orey considered the same small random perturbation for the 
dynamical system with reflecting boundary in smooth domain, and obtained the Freidlin- 
Wentzell's large deviation estimates in C([0, T]; D{A)). They assumed that the coefficients 
are bounded and Lipschitz continuous, the diffusion coefficient is non-degenerate. Using 
the contraction principle, Cepa [5] only considered the large deviation principle of one 
dimensional case based on an explicit construction of the solution (cf. [13]). Multi- 
dimensional case is still open. Compared with the usual SDE, i.e., A = 0, most of the 
difficulties come from the presence of finite variation process K{t). One only knows 
that t I— i> K{t) is continuous, and could not prove any further regularity such as Holder 
continuity. Therefore, the classical method of time discretarized method is almost not 
applicable (cf. [7]). 

Our method is based on recently well developed weak convergence approach due to 
Dupuis and Ellis [S] (see also [21 E])- This method has been proved to be very effective 
for various systems (cf. [HI [TSl IH [ISl [IDl [121 etc.]). Since we cannot prove the following 
uniform estimate as in [TT]: for any p ^ 2 and s, t G [0, T], x, ?/ G D{A) 

sup ¥.\X'{t,x)-X'{s,y)\^'P ^ C{\t- s\P + \x-y\^P), 

^£(0,1) 

one cannot obtain the tightness of the laws of X'^{t,x) in S. Some technical difficulties 
for verifying the conditions (LD)i and (LD)2 below need to be overcome. 

In Section 2, we recall some well known facts about the MSDE and a criterion for 
Laplace principle. In Section 3, we present our main result and give a detailed proof. 
Throughout the paper, C with or without indexes will denote different constants (de- 
pending on the indexes) whose values are not important. 

2. Preliminaries 

We first give some notions and notations about multivalued operators. Let 2^™ be the 
set of all subsets of M"*. A map A : —>■ 2^"^ is called a multivalued operator. Given 
such a multivalued operator A, define: 

D{A) := {x G : A{x) ^ 0}, 
lm{A) := U^^D(A)Aix), 

Gt{A) := {(x, y) G R^'" : X G M*", y G A{x)}. 
We recall the following definitions. 
Definition 2.1. (1) A multivalued operator A is called monotone if 

{yi - y2, xi - X2)js^rn ^ 0, V(xi, yi), (x2, 1/2) e Gr{A). 

(2) A monotone operator A is called maximal monotone if for each {x,y) G Gt{A), 

{y-y',x-x')j^m^O, y{x',y') e Gi{A). 
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Examples 2.2. Suppose that O is a closed convex subset ofW^, and Iq is the indicator 
function of O, i.e, 

^o{x) .- I ^f x^O. 

The subdifferential of la is given by 

dioix) := {yeR"':{y,x-z)^^^0,yzeO} 

0, zf O, 
{0}, if xe Int(O), 
A^, if xe do, 

where lnt{0) is the interior of O and is the exterior normal cone at x. One can check 
that dio is a multivalued maximal monotone operator in the sense of Definition \2. li 

We now give the precise definition of tlie solution to Eq.([T]). 

Definition 2.3. A pair of continuous and (J-'t)- adapted processes {X,K) is called a so- 
lution of Eq. (1) if 

(i) X(0) = X, and for all t ^ 0, X{t) e D{A) a.s.; 
(a) K{0) = a.s. and K is of finite variation; 

(Hi) dX{t) = b{X{t))dt + a{X{t))dW{t) - dK{t), ^ t < oo, a.s.; 
(iv) for any continuous and (J-'t)- adapted processes {a,f3) with 

G Gr(A), VtG[0,+oo), 

the measure 

{X{t) -a{t),dK{t) - P{t)dt) ^0 a.s.. 

We now recall an abstract criterion for Laplace principle, which is equivalent to the 
large deviation principle (cf. [31 HI [15]). It is well known that there exists a Hilbert space 
so that P C V is Hilbert-Schmidt with embedding operator J and {W''{t),k G N} is a 
Brownian motion with values in U, whose covariance operator is given hj Q = Jo J*. For 
example, one can take U as the completion of with respect to the norm generated by 
scalar product 



For a Polish space B, we denote by B{M) the Borel cr-field, and by Ct(B) the continuous 
function space from [0, T] to B, which is endowed with the uniform distance so that Ct{ 
is still a Polish space. Define 

4 := |/i = ^ h{s)ds : h G L\0, T; l^) \ (3) 

with the norm 

where the dot denotes the generalized derivative. Let /i be the law of the Brownian motion 
W in Ct(U). Then 

(4,CT(u),/i) 

forms an abstract Wiener space. 
For T,N >0, set 

Vn ■■= {heel - ^ 
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and 

\ h : [0,T] l'^ is a continuous and (jFi)-adapted 
= I process, and for almost all £ P. ^ ' 

We equip Vn with the weak convergence topology in Then 

Pat is metrizable as a compact Pohsh space. (5) 

Let S be a Pohsh space. A function / : S ^ [0, oo] is given. 

Definition 2.4. The function I is called a rate function if for every a < oo, the set 
{/ G S : /(/) ^ a} is compact in §. 

Let {Z^ : Ct(U) ^ S, e G (0, 1)} be a family of measurable mappings. Assume that 
there is a measurable map Zq : hh^ § such that 

(LD)i For any > 0, if a family {h^,e G (0,1)} C Aj^ (as random variables in Pjv) 
converges in distribution to h E A]^, then for some subsequence e^, Z'''' ' - ' ^ 



converges in distribution to ZqIH) in §. 

(LD)2 For any A^ > 0, if {hn,n G N} C Vn weakly converges to /i G £'^, then for some 
subsequence hn^, Zo^hnJ converges to Zq^H) in §. 

For each / G S, define 

1 

2 {/ie4:7=2oW} 



Hf)--=7,^. ,Jnf_^J\h\\l, (6) 



where inf = oo by convention. Then under (LD)2, /(/) is a rate function. 
We recall the following result due to [3] (see also [TTl Theorem 4.4]). 

Theorem 2.5. Under (LD)i and (LD)2, {Z'^,e G (0,1)} satisfies the Laplace principle 
with the rate function I{f) given by More precisely, for each real bounded continuous 
function g on E>: 



lime log ( exp 



inf{^7(/) + /(/)}. (7) 

/6s 



In particular, the family of {Z'^,e G (0,1)} satisfies the large deviation principle in 
(S,B(E>)) with the rate function I{f). More precisely, let be the law of Z^ in (§, i3(§)), 
then for any B G i3(S) 

— inf /(/) ^ liminf elogz^e(-B) ^ limsupelogz/e(i?) ^ — inf /(/), 

/eB° e^o /e-B 

where the closure and the interior are taken in §, and I{f) is defined by 

3. Main Result and Proof 

We assume that 

(HI) A is a maximal monotone operator with non-empty interior, i.e., lnt{D{A)) ^ 0; 
(H2) cr and b are continuous functions and satisfy that for some Co-, Cb > and all 
x, 1/ G 

\\a{x) - a{y)\\L^^i2.^m) ^ C^\x - y\, 
{x-y,b{x) -b{y))j^„, ^ C^jx - 
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where L2(/^;M"*) denotes the Hilbert-Schmidt space and | ■ | denotes the norm in 



and for some > and n eN 



\bix)\^aii + \xn. 

It is well known that under (HI) and (H2), there exists a unique solution (X^, K"^) to 
Eq.Q in the sense of Definition 12.31 (cf. [E]). Our main result is stated as follows: 

Theorem 3.1. Assume that (HI) and (H2) hold. Then the family of {X'^{t,x), e G 
(0, 1)} satisfies the large deviation principle in § := C([0,T] x D{A)] D{A)) with the rate 
function given by 

where X'^(t,x) solves the following equation: 

dX\t) e biX\t))dt + aiX\t))h{t)dt - A{X^{t))dt, X^(0) = x. 

Remark 3.2. Let D be the half plane in M™, i.e., 

D:={x= {x\x^, x"^) eW^-.x^^ 0}, dD := {x : x^ = 0}. 

Let OId be the subdifferential of the indicator function Id (see Example \2.2\l . Consider 
the following small perturbation of SDEs with refection boundary studied in [1] : 

dX\t) e b{X'{t))dt + ^ea{X'{t))dW{t) - dlD{X'{t))dt, 
X'{0) = xeD, eG(0,l). 

It clearly falls into our formulations and Theorem \3.1\ can be applied to obtain the same 
result as in^. 



For proving this result, by Theorem 12.51 the main task is to verify (LD)i and (LD)2 
with 

§ := C([0, T] X Dp); ^Dp)), = X% Zo{h) = X^. 
This will be done in Lemmas 13.81 and 13.91 below. 

Let G A]^ converge almost surely to /i G Ajj- as random variables in i"^, and 
{X'^'^% K'^'^") solve the following control equation: 

X^'^^(t) = x + [ b{X''''^{s))ds+ [ cx{X''^^{s))k{s)ds 
Jo Jo 

/ cr{X''^^{s))dW{s) - K''^^{t), (9) 
Jo 

which can be solved by Girsanov's theorem, and (X'', K^) solve the following deterministic 
equation: 

X\t) = x+ [ b{X\s))ds+ [ a{X\s))h{s)ds-K\t). (10) 



Let \K\^ denote the total variation of K on [s,t]. We recall the following result due to 
Cepa [6] (see also [HI Propositions 3.3 and 3.4]). 

Proposition 3.3. Under (HI), there exist a G ]R'",7 > 0,yU ^ such that for any pair 
of (X, K) with the property (iv) of Definition \2.3\ and all ^ s < t ^ T 

(X(r) - a, dK{r))^^ ^ ^\K\t - /i T |X(r) - a\dr - 7/^(t - s). (11) 



Moreover, for any pairs of {X,K) and {X,K) with the property (iv) of Definition \2.3i 

(X(t) - dK{t) - dK{t))^^ ^ 0. (12) 
Using this property, we first prove the following uniform estimates. 
Lemma 3.4. For any p ^ 1, there exists Cp^T,N > such that for any e G (0, 1) and 

x,ye D{A) 



El sup \X^''^it,x)-X^'''^it,y)\'^] ^C,,T,N\^-y\''P. 



|2p 



(13) 



Proof. Set 
and 



Z,{t) ■.= X^''^^{t,x)-X^''^^{t,y) 



Ais) ■.= a{X^''^^is,x))-criX^''^^is,y)). 
By Ito's formula, (H2) and ( fT2l) . we have for any p ^ 1 



+2p / \Z,{s)\'P-'{Z,{s),A{s)K{s))^^ds 
Jo 

+2p^e f \Z,{s)\''^-\Z,{s),K{s)dW{s))^^ 



-2p J \Z,is)\'^-'{Z,{s),dK^''^{s,x)-dK^'''-is,y))^^ 
-pf\Us)\-- (llA(.)ir + 2ip - ,) (^^(^)^A(.)A^(.)Z.(.)), 



ds 



^ |Z,(0)|2P + C/ \Z,{s)\^Pds + 2p \Z,{s)\'P ■ \\K{s)\\i2ds 
Jo Jo 

+2p^e f \Z,{s)\'^~^Z,{s),K{s)dW{s))^^. 



Set 

^7(t):=E( sup |Z,(s)|2p 
\se[o,t] 

By BDG's inequality and Young's inequality, we have for any 5 > 



E 



sup 2p^e / \Z,{s)\^^-\Z,{s),A{s)dW{s))^^ 
t'e[o,t] Jo 

^CE (£|^e(^)r^-1|A(^)*^e(^)||?.d.) 

^ CE I sup |^,(s)pP f |Z,(s)pPds 
Yse[o,t] Jo 

^6-g{t) + Cs fE\Z,{s)\^Pds. 
Jo 



1/2 



(14) 



Similarly, we have 



E 



2p / \Z,is)\'^ . \\kis)h.ds 



^ CNE (^j^ \Z,{s)\^Pds^ 

^ S-git) + Cs,N f E\Z,{s)\^Pds. 
Jo 



(15) 



Letting 5 = 1/4 in f|T4|) and f|T5|) and combining the above calculations, we get 

9it)^9{0) + l:g{t) + C fE\Z,{s)\''^ds^2g{0) + 2C f g{s)ds, 
^ Jo Jo 

which gives the desired estimate by Gronwall's inequality. 



□ 



Lemma 3.5. For any p ^ 1 and x G D{A), there exists Cp^T,N,x > such that for any 
ee (0,1) 



E sup \X^'''^{t,x)\^P] +E\K'''''{-,X)\^T ^Cp,T,N,.. 
\te[o,T] J 

Proof. First of all, as in the proof of Lemma [3.41 we may prove that 

e( sup |X^'"^(t,x)pM ^C,,T,N,.. 
V*e[o,T] / 



(16) 



(17) 



Let a G lnt{D{A)) be as in Proposition 13.31 By Ito's formula, flTT]) and (H2), we have 
hx''''^{t)-a\^ = i|x-a|2+ / (X^''^^(s) -a,6(X^'^^(s)))jj™ds 



Jo 

+Ve [ {X^'''{s)-a,a{X^''^^{s))dW{s))^ 
Jo 



(X^'^^(s) -a,dir^''^^(s)) 



+ '-j\w{x^Hs) 



<: -|x-a|' + a/ \X''^^{s) - a\'ds 



+ / (X^''^^(s) -a,6(a))j5„ds 
Jo 

+N (^j^ |(t(X^''^^(s))*(X^'^^(s) - a)\Ms^ 



1/2 



/ {X^'^'^is) - a, a{X''''^{s))dW{u))^^ + fi-ft 



Jo 

+^ j\\x'Hs)\+^m'ds. 
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The desired estimate now follows by (fT7|) . 
Define 

w,{t,x) := [ a{X^{s,x)){h,{s) -h{s))ds. 
Jo 

The following lemma is easy by Ascoli-Arzela's lemma. 



Lemma 3.6. We{-,x) converges a.s. to zero in C{[0,T], D(A)). 

We now prove the following key lemma. 
Lemma 3.7. X^'^" defined by Ij^ converges in probability to defined by [W\) 
Proof. Set v,{t) ■=v,{t,x) :=X^'^^(t,x) -X^{t,x). Then 



v,{t) = K''''^{t)-K\t)+ / {b{X''''^{s))-b{X\s)))ds 







+ / ia{X''''^{s))k{s)-a{X\s))h{s))ds 
Jo 

+ / (T(X^'^'(s))diy(s). 

Jo 



By Ito's formula, we have 



+2 / {v,is),biX^''^is))-biX'is)))^^ds 



+2 / {v,{s),{a{X^''^{s)) - a{X\s)))h,{s))^^ds 







+2 / {v,{s),aiX\s))iKis) ~ h{s)))^^ds 



+2v^ / {t 
Jo 



v,{s),aiX^''^^is))dWis)\ 



+e / MX^'''{sm',.^^^,ds 







=: m + m + mt) + m + m + m. 

It is clear that by <^ 

im ^ 

and 

m^^ct I \v,{s)\^ds. 

Jo 

By BDG's inequality and (H2) we also have 

E I sup \Il{t)\\ + E I sup |J^(t)| I ^ C^e. 

W[0,T] / \t(i[0,T\ / 



As estimating (fT5|) we have 



E ( sup |/|(s)| ) ^ ( sup \v,{s)\A + C [ E\v,{s)\Ms. 
\s6[o,i] J 2 \^se[o,t] y Jo 



Set 



g{t) :=E I sup 

se[o,t] 



Then we have 



^7W^^^?W + C^v^+e( sup \Il{s)\] + C [ E\vM\Ms, 
^ \se[o,t] j Jo 



which imphes that 



g{t) ^CV^ + 2eI sup \Il{s)\]+C [ g{s)ds. 

\se[o,T] J Jo 

By Gronwall's inequahty we get 

e( sup \v,{s)\A ^Cv^+Ce( sup 1/4(5)1). (19) 
V^e[o,T] J \se[o,T] y 

We now deal with the hard term II. By Ito's formula again, we have 







'w,is){b{X^'^^{s))-b{X\s)))ds 







- I w,{s){a{X''^^{s))k{s) - a{X\s))h{s))ds 







-v^ j w,{s)a{X''''^{s))dW{s) 
For any 5 > and i? > 0, we have 



any 5 > and i? > 0, we have 

P I sup \Il^{t)\ ^5] = P I sup |/4\(t)| ^ 5; sup |?;,(t)| < ) 
ye[o,r] y \te[o,T] ie[o,T] y 

-^P I Riir, I /'Iff ^ I > A- Riin U) r/^l > /? 



+P ( sup \Il^{t)\ ^ 6; sup \v,{t)\ ^ R 
^te[o,T] te[o,T] 

^ P I sup \w,{t)\ ^ 6 / r] + P { sup \v,{t)\^R 



By Lemma 13.61 and (1161) we know 



limP sup \Ili{t)\ ^5=0. 



^ie[0,T] 
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Noting that 

sup sup \w,{s)\-{\K^''^{-)\'t+\K'{-)\'^), 

te[0,T] s€[0,T] 

as above, we also have 

sup 1/42(^)1 ^ in probabihty. 

te[o,T] 

Similarly, we have 

sup 1/43(^)1+ sup 1 744 (t)| ^ in probability. 

te[o,T] te[o,T] 



Moreover, by BDG's inequality we have 

sup \l 



V^E sup \Il,{t)\ I ^ CV^. 
\te[o,T] 



Combining the above calculations, we get 



sup 1/4(^)1 ^ in probability. 

ie[o,T] 



It is easy to see by ( 1161) that 



Hence 



Substituting this into (HM we obtain 



sup E sup \Il{t)\'^ I < +00. 

ee(0,l) \t6[0,T] 



limE I sup \Il{t)\ I = 0. 

te[o,T] 



limE sup \v,{s)\^ = 0. 



Thus, we have proven that for all x E D{A) 

sup — s> 0, in probability. 

i6[0,T] 

We now strengthen it by Lemma [3.41 to 

$,n,e '■= sup — s> 0, in probability. 

te[0,T],xeD{A),\x\^n 

Set 

:= D{A) n{xeW : \x\ ^ n} n 511^, 

where 5 > and 5Z™ denotes the grid in with edge length 5. It is clear that there 
are only finite many points in D^. Hence 

in^:= sup ^ 0, in probability. 

ie[o,T],xeD« 

For any x G M™, let xs denote the left-lower corner point in 517^ so that 
Noting that 

^ 2^^_^ + 2 sup sup \v^{t,x) - V^{t,Xs)\^, 

xeZ)(A),|xKn*e[0,T] 

10 



we have for any (3 > and some a > 

P{^n,e > 4/3) ^ P [C,e + SUp \v,{t,x) - V,{t,Xs)\^ > 2(3 

\ xeD(A),|x|<nte[0,T] 

\xeD(A),|x|<n*e[0,T] 

^ P{^1,>(3)+eI ^p sup \v,{t,x)-v,{t,xs)n /P 

\x&D{A),\x\<inte[0,T] J 

where the last step is due to Lemma 13.41 and Kolmogorov's criterion. First letting 6 be 
small enough, then e go to zero, we then obtain 

limP(en,. >4/?) = 0. 

which yields the desired result. □ 
Lemma 3.8. (LD)i holds. 

Proof. Let be a sequence in Aj^ converge to h in distribution. Since Pat is compact 
and the law of W is tight, {h^, W} is tight in P^r x Ct(U) by the definition of tightness. 
Without loss of generality, we assume the law of {he, W} weakly converges to /i. Then 
the law of h is just /i(-,Cr(U)). Indeed, for any bounded continuous function g on V^, 
we have 

E{g{h)) = lim E{g{he)) = [ (?(/i)/i(d/i, Ct(U)). 

By Skorohod's representation theorem, there are {Cl,P) and {h^jW^} and {h,W} such 
that 

(1) {he,W^) a.s. converges to {h,W); 

(2) {he, W") has the same law as {h^, W); 

(3) The law of {h, W} is /x, and the law of h is the same as h. 
Using Lemma [37ft we get 

$ ( ^ / heds + wA ^ X'', in probability. 







where $ is the strong solution functional (cf. ^). From this, we derive 



$ ( ^ / heds + W] ^ X'', in distribution. 



Ve JO 

Thus, (LD)i holds. □ 

Similar to the proof of Lemma 13.71 one can easily verify that 
Lemma 3.9. (LD)2 holds. 

Thus, by Lemmas 13.81 13.91 and Theorem 12.51 we have proved Theorem 13.11 
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